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$h ' Abstract 

We consider higher-derivative perturbations of quantum gravity and quantum field theories in curved 
space and investigate tools to calculate counterterms and short-distance expansions of Feynman diagrams. 
In the case of single higher-derivative insertions we derive a closed formula that relates the perturbed 
one-loop counterterms to the unperturbed Schwinger-DeWitt coefficients. In the more general case, we 
classify the contributions to the short-distance expansion and outline a number of simplification methods. 
Certain difficulties of the common differential technique in the presence of higher-derivative perturbations 
are avoided by a systematic use of the Campbell-Baker-Hausdorff formula, which in some cases reduces 
the computational effort considerably. 
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1 Introduction 

In quantum gravity infinitely many independent couplings are necessary to remove the diver- 
gences. Practical tools to do systematic computations with non-polynomial lagrangians are not 
available today. In this paper we investigate techniques to express the one-loop counterterms of 
the most general lagrangian in closed form. Consider a generic gravitational action constructed 
with the curvature tensors and their derivatives, 



The one- loop counterterms are collected in a functional Sidiv uniquely determined by Stree, 



If S'ldiv (5'tree) could be written explicitly, it would be possible to search for special Stree 's con- 
taining a finite number of independent couplings, such that 



up to two-loop corrections, where S[ Tee coincides with Stree up to redefinitions of fields and 
couplings. If Stree were so special to satisfy analogous identities for the two-loop and higher- 
order counterterms, then it would define a "renormalizable" theory. 

The search for renormalizable theories beyond power counting is not an easy task, but can 
teach us a lot about the structure of counterterms and their classification. In quantum gravity 
renormalization turns on vertices with dimensionalities greater than four. It is well-known that 
in the absence of matter, the one- loop counterterms can be eliminated with a field redefinition of 
the metric tensor [1] . The first new vertex is cubic in the Riemann tensor and removes a two- loop 
divergence [2j [3] . The corrected quantum gravity lagrangian reads 



Expanding (jl.ip around a background metric, the one-loop Feynman diagrams are encoded in 
the determinant of a differential operator containing higher-derivative terms. In general, the 
higher-derivative terms can be treated perturbatively or non-perturbatively. In the former ap- 
proach [HO [3] ("quantum gravity") they are viewed as perturbations of the Einstein lagrangian: 
the theory is non-renormalizable, but perturbatively unitary. In the latter approach [H EJ [6] 
("higher-derivative gravity") they are used to improve the behavior of Green functions at short 
distances: the theory is renormalizable, but not unitary. Here we are interested in the former 
approach, which is equivalent to study the insertions of higher-derivative operators in the Feyn- 
man diagrams of quantum gravity. Observe that in (jl.ip higher powers of the curvature tensor 
generate perturbations with an arbitrary number of derivatives. 
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We illustrate our techniques in the case of a scalar operator of the form 



oo 



H = H + H 1 



(1.2) 



n=0 



For our purposes H\ can be treated perturbatively. We loose no generality if we assume that the 
tensors V^ 1 '"^" are completely symmetric. Indeed, commuting the covariant derivatives every 
antisymmetric component of V^ 1 '"^ 1 can be reduced to a combination of V-terms with fewer 
indices. We investigate tools to study the perturbative expansion of the .ff-determinant and 
simplify the computation of its coefficients. Our arguments are general and their extension to 
spinors, spin-1 fields and the graviton is direct. The case of gravity is addressed. 

Calculations in quantum gravity are conveniently done using the background field method 
[8] and, at the one-loop level, the Schwinger-DeWitt techniques [H1CE0], because they manifestly 
preserve covariance. The common approach is to perform a Schwinger-DeWitt expansion of the 
Green function, derive a differential equation for its coefficients and work out their short-distance 
expansion by repeated differentiation. However, in the presence of higher-derivative perturbations 
the differential approach has some difficulties, which can be easily overcome working at the level of 
operators and systematically using the Campbell-Baker-Hausdorff (CBH) formula before taking 
the coincidence limits. In some cases this approach simplifies the calculation considerably and 
allows the derivation of some closed formulas. Moreover, it singles out that a number of involved 
expressions are just total derivatives and so can be neglected for the purposes of renormalization. 

The paper is organized as follows. In section 2 we recall the Schwinger-DeWitt approach and 
explain the difficulties of the differential technique. In section 3 we introduce the CBH approach 
and derive the closed counterterm formulas (|3.10p and f|3. 1 1 [) for the most general single insertions. 
We also prove the new identity (|3.9p . In section 4 we describe the method in general and make 
explicit computations. In particular, formula (|4.10p for three-derivative perturbations is a new 
result. In section 5 we describe how the techniques apply to gravity. Section 6 contains our 
conclusions. 

2 Difficulties of the differential approach 

Given an operator H, define the function H(s;x,x'), s > 0, as the solution of the equation 





with the boundary condition 
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If H has for example the form 

H = □ - £R, (2.3) 

where □ denotes the covariant D'Alembertian, the Schwinger-DeWitt expansion of the associated 
function Ho(s; x, x'\ £) reads 

tf o 0w';O = "74^ ex P (^|^) f^{is) n A n {x , x 1 ^) , (2.4) 
with the boundary condition 



lim A (x,x';£) = 1. (2.5) 

a:'— he 

In (|2.4p <t(x, x') is one half the squared geodesic distance between x and x' and satisfies 

cj, a(x,x) = a ]lx {x,x) = 0, cr ;tiu (x,x) = g^ u (x). (2.6) 



2 ,M 



Equation (12. ip . with — > i7o> generates a differential recursion relation for the coefficients A n , 
n > 0, namely 

(n - 2)A n + a'^An^ + -A n U x a = (D x - £R) A n _i, (2.7) 



with A^i = 0. By repeated differentiation, the recursion relation (|2.7p can be used to calculate the 
short-distance expansion of the Schwinger-DeWitt coefficients A n (x, x'; £). For this purpose, it is 
sufficient to compute the coincidence limits A n (x,x;^), which are called "diagonal coefficients", 
and the coincidence limits of the covariant derivatives of A n (x, x'; £), which are called "off-diagonal 
coefficients". The calculational method just described will be called the "DeWitt differential 
approach". The coincidence limits will be denoted with an overline. The first two coefficients, 
Ai and A2 have been computed by DeWitt in [10j . A3 by Sakai in pTJ and Gilkey in [12], A4 by 
Amsterdamski, Berkin and O'Connors in [13] and Avramidi in [14j, A^ by van de Ven in [15]. In 
[16] the first off-diagonal coefficients have been recently worked out to a considerable order and 
number of derivatives. 

The one-loop contributions to the generating functional V of one-particle irreducible functions 
read 

r (1) = --y 5 -J & A x^(x)H{s-x,x-t). (2.8) 

Although our techniques are general, we focus on the scheme-independent (logarithmic) diver- 
gences. In the notation commonly used in dimensional regularization (In 5 — > — l/(2e)), we have 



Thus, to study the one-loop renormalization one has to calculate the coincidence limit of 
the second Schwinger-DeWitt coefficient. In the DeWitt differential approach, this goal can be 
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achieved repeatedly differentiating equation f|2.Tj) and the first of (|2.6p . and taking coincidence 
limits with the help of (f23|) and (pUj) . 

However, the differential approach is not convenient to study higher-derivative perturbations. 
The reason can be appreciated already in flat space. Consider 

H = H + Hi, H = d 2 , H l = \{d 2 ) 2 . (2.10) 

The unperturbed flat-space Green function reads 

Ho{s; X ~ X) = ~ (4^ 6XP J • 

The first observation is that the Schwinger-DeWitt expansion (|2.4j) of H (s; x — a/) needs to be 
replaced with a sum containing arbitrary negative powers of s, namely 

E M n Bn(Kz). (2.11) 

Nevertheless, at each order in A the sum is bounded from below. In particular, at 0(A) the sum 
starts at n = —3. To this order, equation (|2.ip gives the relations 



Differentiating these relations a suitable number of times and taking the coincidence limits (x 
0), we find B_2 = -B-3 = O-B-3 = 0, plus the relations 



a^B^ 3 = 12A, B-i + DS_2 = -6A, (2.12) 

which are valid up to higher orders in A. Two equations give the first of (|2.12p . so one quantity, 
□1?_2, remains undetermined. 

More generally, the recurrence relations for the coefficients B_k with k > have the form 

kB_ k -x^d^B_ k = P k + 0(X), (2.13) 

where P k possibily depends on d k '~ k B- k ' with k' > k. The left-hand side of (|2.13p vanishes, in 
the coincidence limit, when k derivatives act on it. Therefore (|2.13p does not provide information 
about d k B_ k . This ambiguity has the following explanation. The initial condition (|2.2p deter- 
mines the solution uniquely. While in the unperturbed problem (|2.2p is exhaustively expressed 
by (|2.5|) . in the perturbed problem it is expressed by Bq = 1 plus suitable relations among the 
BnS with n < 0. The extra relations, however, are not immediately readable from the expansion 
(|2.1ip . because of the negative powers of s contained in the sum, and need to be worked out 
independently. 
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For these reasons it is more convenient to pursue a strategy that incorporates the boundary 
condition (12. 2p automatically. This goal is achieved writing 

H(s;x,x') = (x\e iSs \x'), (2.14) 

where | x) are position eigenstates, x^ \ x) = x^ \ x), {x' \ x) = 5(x' — x). Noting that in the case 
(|2,10p Hq and H\ commute, we can write 

H(s;x,x') = (x | e iS ° s e iS ' s | x') = ^ s H Q {s;x,x') = --^—j^d 2 ?^ ( ^ ~ x') 2 \ 
This procedure does not contain any ambiguity, and easily leads to 



B_ x = 6A, UB-2 = -12A. 

In the rest of the paper we use this strategy in curved space. We name it "CBH approach" , 
because it involves a systematic use of the CBH formula. Besides avoiding the difficulty just men- 
tioned, in some cases the CBH approach reduces the calculational effort considerably. Moreover, 
it allows us to calculate each coefficient B n directly, without having first to recursively calculate 
the U m 's with m < n. 

3 The CBH approach 

In curved space, formulas (12. 14ft and (|2.1ip are replaced by 

H(s;x,x';Z) = (-g(x))- l ^(x \ ^ | x'>(- 5 (x'))- 1/4 
i ( ia(x, x' 



(4tt) 2 s 2 



exp(^g^) f; (isrB n (x,x';0, (3.1) 



where H = (- 5 ) 1 / 4 J §(-c/)- 1 / 4 . Define H and Hi as in JL2|) and write H = H Q + Hi. The CBH 
formula reads 

jHs = e iH s J2 S^L / dCl • • • dCn T [^(Cl) • • • 5n(Cn) , (3.2) 

where T denotes the ordered product and 

Fx(C) = e-^H^ ^ = Y (~ ls p n (adg ) w gi- (3.3) 

z — ' n! 

n=0 

with (ad^4)-B = [A, B]. Consider, for example, the first order in Hi, namely the diagrams that 
contain a single insertion of the perturbation. We have 

H{s;x,x';Z)=H {s;x,x';Z) + is f d({-g(x))- 1/4 (x | e ^(i-0# ie ^o< | x '){- g tf))-V* 



o 



+0(H 
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The one-loop contributions (|2.8p to the T functional become 



= - / ds d( d 4 x (x I e iH '(l-C)H ie iH , 



X). 



The ^-integrand is just the trace of the operator contained between the bra and the ket. We can 
use the ciclicity of the trace and get 



r« = - / ds / d 4 x (x | H ie 



JHqs 



ds 



HiH Q {s;x,x';£) 



(3.4) 

Thus, to compute the one-insertion one-loop diagrams it is sufficient to act with H\ on the 
unperturbed function Hq and then take the coincidence limit. The divergent part is given by the 
0(1/ s) contributions to the square bracket in (|3.4p . namely 



r 



(i) 

div 



4e 



d A xy/-g(x) HiHo(s]x,x']t) 



(3.5) 



where the superscript s _1 is to emphasize that only the coefficient of 1/s has to be kept, after 
inserting the Schwinger-DeWitt expansion for the unperturbed function Hq(s; x , x' ; £) . 

To illustrate these facts in a simple example, consider a complex scalar field ip in curved space, 
described by the lagrangian 

-£= = -d^gTduip - Zmpp + ¥ (V + VD^ + V^D^Dv + V^ P D„D V D P + •••)¥>, (3.6) 



where all tensors V pv "' are symmetric. 

The one-insertion divergent terms are then 



1 



r dii = ^ / d 4 *v^ [{V + V^D^ + V^D^D U + V^D^Dp + • • •) H (s; x, x'; £)] " x , =x 



(3.7) 



The first two types of terms give immediately 



div 64tt 2 £ 
The two-derivative term gives 



d 4 x^ (V^Ai + 



(1 ~ 6g) 
24(4vr) 2 e 



.(1) 
div 



1 



64vr 2 e 



d 4 xy^ ( - -V£A 2 



960( 4 7r )2 £ / d^v 7 ^ (DiV + (1 - 100#V + (3 - 20£)i^ - 2R^R pii(TV ) 

-V* (2R vp R vp + 4(1 - 5£)OR + (60£ 2 - 20£ + l)i? 2 )] . 
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The three-derivative term gives 



1 



647r 2 e 



& A XyJ—g 



— I R ^ V up + — ff w /? V" ra -I- — R^ v R V a/3 
4q ] >m ^ p 20 V/» k k> 80 v ol;P 



1 e , 3, 2 



80-1 + 4^ 



(3.8) 



The four-derivative term gives 



The coincidence limits A3, A^u and A\-^ upcf that are necessary to write this expression explicitly 
have been worked out in [16] and rederived by ourselves. 
The five- and six-derivative term gives 

respectively. Apart from A4 , which has been computed in [13] and [2], the unperturbed coeffi- 
cients appearing in these formulas have not been written in the literature. 

A useful identity. The formula for can be simplified using the identity 

^a^i-mu) = °' Vn>l, (3.9) 

where the parenthesis means complete symmetrization. Expressions such as V^'^a-^...^ and 
similar are thus identically zero. This property reduces the number of cr-derivatives that need to 
be computed to work out T$ . Formula (|3.9[) can also be used to derive the coincidence limits 
cT;^...^ in a more efficient way. 

The proof of f|3.9|) can be done by induction. For n = 2 the identity is true, since <7 ;/il/i2At3 = 0. 
Assume that it is true up to n = n > 2. Taking one derivative of the first equation of (|2.6p . we 
get 

Now, take n + 1 derivatives of this equation and symmetrize completely in those. We get 

k=0 \ K ) 
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where A and \i are excluded from the symmetrization. Now, take the coincidence limit of this 
expression and use the inductive hypothesis, together with ct^J = 0. The result simplifies to 



C-//Y//.1 ,1 =(n + l)cr' A , a.,.\ + a' x , ,a : , 



Using o^7 = fif//i/ we arrive immediately at 



which proves the statement. We have checked (|3.9p explicitly up to n = 7 included, using the 
complete expressions of cr ;Atl ... /Jn , n < 8, derived with a computer program. 

The general formula. Using (13. 9p it is possible to work out the general formula 

oo [n/2] 



F - = 64A / d4 ^E E ^Kn -»)l ^-^- trtV '°""-- <3 - 10) 

n=0 k=0 V ' 



where [n/2] is the integral part of n/2 and tr^U means that k pairs of V-indices are traced. The 
formula is derived as follows. Consider (|3.7|) with the perturbation V J,1 '" IMl D lxl ■ ■ ■ D^ n : 

r£ = 4^ / d^V^lv^D^.-.D^ exp (^f^) f>)»- 2 ^(*y ; fl 

Since the derivatives are symmetrized, any time three or more of them act on a(x, x') the contri- 
bution vanishes in the coincidence limit. Moreover, since a-n vanishes, only two derivatives can 
act on the same a, all others having to act on the A n 's. Since = g^, two derivatives acting 
on a trace a pair of U-indices. The remaining combinatorics are then straightforward and give 

_i 

J ^ fc=0 V ; m=0 

which proves (|3.7p . 

In renormalization theory a further simplification applies. Indeed, it is not necessary to include 
in (|3.6p independent terms proportional to □</?, because they can be converted into terms of other 
types by means of c/9-field redefinitions. Up to terms proportional to the <£>-field equations, □</? 
can be replaced with 

SUp -(y + V^D^ + V^D^Dv + V^D^D u D p + • • >. 

Thus, a repeated use of field redefinitions can eliminate all boxes acting on <p. Moreover, any 
couple of contracted U-indices can be moved to a box acting on <p commuting the covariant 
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derivatives, up to V-terms with fewer indices. Thus, it is sufficient to take symmetric, traceless 
Vs. The final result is then just 

. oo 

T * = 64^ / d 4 ^E^^: V*"^. (3.H) 

J n=0 

4 Calculation of total derivatives and multiple insertions 

In a variety of computations, for example the short-distance expansion of Green functions, to- 
tal derivatives have to be kept. Moreover, with multiple -Hi-insertions, even neglecting total 
derivatives, the calculation simplifies much less than with single i^i-insertions. In this section we 
describe how the calculation based on the CBH approach proceeds in the general case and report 
a number of sample and new calculations of the perturbed coefficient Bi- First we focus on the 
total-derivative corrections to the single i?i-insertion results computed in the previous section. 
Later we classify the structure of contributions in the general case. 

The scalar-potential, one-derivative and two-derivative results (|4.4|) . (|4.5|) and (|4.7|) are known. 
They can be derived in a variety of conventional ways. We rederive them with our techniques to 
illustrate the CBH approach. They are also useful to introduce the more difficult derivation of 
the three-derivative new result (|4.10p . 

Write 

H(s; x, x'; = H (s; x, x'; £) + H^s; x, x'; £) + O(Hf). 
From the CBH formula (|3.2|) we get 

ffifos.a'ifl = V [ ^—r ( a dH ) n HiHo(s-,x,x'-,0- (4-1) 
^ (« + l)! 

Suppose that the interaction H\ contains at most m derivatives and that we are interested in the 
calculation of the coincidence limit of a given perturbed coefficient, say B^. Since Hq contains 
at most two derivatives and each commutator ad-ffo raises the number of derivatives by one 
unit, (adHo) n Hi contains at most m + n derivatives acting on the function Hq{s;x,x'). When 
derivatives act on the exponential prefactor F = exp (a(x,x')/(2is)), they lower the s-power. 

Although each derivative acting on F lowers the s-power by one unit, to give a non-trivial 
contribution in the coincidence limit derivatives have to act on F at least in pairs, because = 0. 
Thus, in the coincidence limit, the s-power can be lowered by at most [(m + n)/2] units. The 
lowest s-power that multiplies the unperturbed coefficient A^i is 

+ k' - 2. (4.2) 



n + 1 



m + n 
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Here the factor s n+1 comes from (|4.ip . while s k ~ 2 multiplies A^i inside Hq(s; x, x'). Now, in the 
Schwinger-DeWitt expansion (|3.1|) of the function H(s;x,x';£), the coefficient is multiplied 
by k — 2 powers of s. Equating this number (|4.2p . we see that, for the purpose of computing the 
single-insertion perturbations to B^, the sum in (|4.ip becomes finite. It stops at the n such that 



n + 1 



m + n 



k. 



(4.3) 



With multiple insertions, say j, the sum of (|3,3p is raised to the power j. Call n the total power 
of ad.f/o contributing from the i?i(C)s and rrij the total number of derivatives carried by the j 
insertions. Then equation (14. 3D is generalized to 



j + n 



rrij + n 



fc, 



The list of contributions stops when the total power of adi?o reaches the value n. 

This counting proves that the CBH method is consistent with the perturbative expansion, and 
the coincidence limit of each perturbed Schwinger-DeWitt coefficient can be calculated algorith- 
mically. However, the calculation can become lengthy quite soon, even for computer programs. 
We now compute the single-insertion perturbations to B% for the cases considered in the previous 
section and classify the contributions of multiple insertions in more detail. 



Scalar-potential perturbation. The simplest perturbation is the scalar potential V(x). In 
formula (|4.ip the term with n = has one power of s, so it gives a contributions proportional 
to A\. The term with n = 1, contains two powers of s and at most one derivative acting on 
Hq(s;x,x'), so its contributions are proportional to Aq-^, which vanishes, and Aq = 1. The term 
with n = 2 contains three powers of s and at most two derivatives acting on Hq, which lower 
the s-power by at most one unit when they act on the prefactor F. This contribution is again 
proportional to Aq, and OV. The terms with n > 3 do not contribute, because they contain too 
many powers of s and too few derivatives to lower them. Working out the commutators we get 

A^= Q-e) W + W. (4.4) 

One- derivative perturbation. Now, consider the perturbation V^D^. The first term of (|4.ip 
gives a contribution proportional to Ai-n. The second term contains two powers of s and at 
most two derivatives acting on Hq, thus it gives contributions proportional to ^4o;2 and A%, where 
j4o;2 denotes any object with less than three derivatives, namely Aq-^, Ao ]fl and Aq. The third 
term has at most three derivatives on Hq. Two of them are used to lower one s power, so the 
contribution is proportional to Aq-\. The third term has four derivatives that have to be used 
to lower the s-power by two units, giving a contribution proportional to Aq. The terms of ()4.ip 
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with n > 4 do not contribute, because they have too many s's and too few derivatives acting on 
Hq. The result is 

&B~2 = -^^-RV* - ^UV*. (4.5) 

We know that the terms with n > 1 are total derivatives. Using the differential approach it is 
not easy to recognize the presence of such total derivatives, which are often very involved. The 
simplest of them is 

(x | e jHos (ad^o)-H : i I x). 

Let us inspect it more closely, to see what kind of relations it generates. First observe that we can 
always move the covariant derivatives away from V p , eventually adding other total derivatives. 
When we do this, we get a relation of the form 

V^Jfj, = total derivative, 

for some current J**. Next, integrating this relation over spacetime and using the arbitrariness 
of V^, we obtain the identity = 0. Finally, substituting the ^-coincidence limits, we obtain a 
relation for the ^-coincidence limits. The result is 



= D„At + tA Q R p + nA Q;fl - 2D V A^ V + A . U R». (4.6) 
Notice that some derivatives are taken before the coincidence limits, others are taken after the 



coincidence limits. The values of Aq-^ v , Aq-^, Aq and Ai are reported in the appendix and indeed 
satisfy (|4.6H . More complicated identities are generated by the other terms of (|4.1|) . 

Two- derivative perturbation. Let us consider the perturbation V^ u D^Dy. The term with 
n = in (|4.ip contains two derivatives, that can either act on the exponential prefactor of Hq, 
lowering the s-power by one unit, or on the unperturbed Schwinger-DeWitt coefficients contained 



in the expansion of Hq. The resulting contribution is a linear combination of Ai and A\.^. The 
term with n = 1 contains at most three derivatives acting on Hq, two of which can act on the 



exponential prefactor. The result is a sum of A±-i and Aq-^. The third term of (|4.ip contains at 
most four derivatives on Hq. Four or two of them can lower the s-power by two units or one, 



respectively. The contributions of this term are proportional to A\ and Aq-2- Similarly, the terms 



with n = 3 and n = 4 give contributions proportional ^4o;i and Aq, respectively. 
The final result is given by 



- — —RVZ + - — —DJR^V^) + —R^nV^ + —R V U V.Z + — - 
18 ,fH/ 30 y ^ ' 36 M 90 ' up 15 



30 ^ 60 M 90 ' * 45 M 9 20 

-—RuvR^Vpa ~ ^—^R-utrV^ - — R„^trV ; ^ - ^RDtrV - — D 2 trV 
45 ^ p 60 4 180 M 72 120 
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(4.7) 



m- derivative perturbation. In the general case, namely a perturbation V pl "' pm D^ 1 ■ ■ ■ D^, 
the first contribution of (14. lh gives the list of terms written in (|3.10j) . Each commutator with Hq 
in (|4.ip raises the s-power by one unit and the number of derivatives by one unit. If the new 
derivative does not act on the exponential prefactor of Hq, we have 

-» (4.8) 

If the derivative acts on the exponential prefactor, then it must absorb a second derivative, to 
give a non-trivial contribution. In this case, both the s-power and the number of derivatives are 
lowered by one unit: 

Ai~ Ak-j-i. (4.9) 
Combining the two operations the contributions fit into the following scheme: 



n 


= 


A[m/2}+l;a{m) 


A[m/2};u(rri)+2 


Ai-m 




n 


= 1 


A[ m /2]+l;cr(m)-l 


A[ m /2};a(m)+l 


Al-m-l 


Ao ; rn+l 


n 


= 2 




A[m/2];a(rn) 






n 


= 3 




A[ m /2];cr(m)-l 






n 


= m 






A X -Q 


A -2 


n = 


m + 1 








A -i 


n = 


m + 2 








A -o 



Here a(m) = if m is even, a{m) = 1 if m is odd. Observe that the coefficient Ao ;m +2, the most 
involved of all, does not contribute. 

For example, for m = 3 we have contributions 



n 


= 


A 2 ;l 


Al;3 




n 


= 1 


A~ 2 


A 1;2 


A -A 


n 


= 2 




A 1;1 


A ;Z 


n 


= 3 




M 


Aq;2 


n 


= 4 






Aq-i 


n 


= 5 









Each of these coefficients are available in the literature and have been recalculated independently 
by us. The m = 3 perturbed coefficient reads: 

AB2~ = — ( -n 2 + -RuvW + -BD + -R 2 + aR + R.^ + -R^D^D U ) V£Z 
40 V 2 2 M 6 4 ,M 3 M J p ' 
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l(-R + -a] VZ P „ ~ —R-uVZ p - —KV-ZZ - —R U uR wca V^ ra - —R, u av.^ a 
g m 5 I 20 60 60 24 



-| (i?D + (1 - 30ii 2 + QR+ ^ t £ M ) ^ + | (RR^ U + fl.^) Kf a . (4.10) 

The unperturbed coefficients necessary for the m = 4 result exist in the literature |16j . For m > 4 
the necessary coefficients can be derived with computer programs, but an increasing amount of 
time is required. 

Squared- Laplacian perturbation. An interesting case is the perturbation An 2 , which is a 
linear combination of the perturbations V^'^D^ ■ ■ ■ D^ m with m = 1, 2, 4. The commutators 
are much simpler in this case and the final result is 



AB 2 = A£ 



—RnR^ + —R^Rau + —a 2 R + — RnR + £ ( I _ £ ) i? 3 
15 ^ 45 ,M 30 18 \ 6 ' 



Multiple insertions. So far, we have classified the single i^i-insertions, but the analysis can 
be generalized to more insertions. We do it for the calculation of a generic perturbed coefficient 
-Bfc. Each multiple- insertion contribution is made by a certain number of Hi's and a certain 
number of (ad-ffo)'s acting on them. Call n the "level" of the contribution, namely the total 
number of HQ-comm\itaX>ois. Denote the total number of f/i-insertions with r. In each insertion, 
pick a perturbation V lll " ^ m D fjLl ■ ■ ■ D^ m , not necessarily with the same number m of derivatives. 
Call d the total number of derivatives carried by such perturbations. Since Hq has two derivatives 
at most, the total number of derivatives acting on Hq is at most equal to d + n. Acting on the 
exponential prefactor of Hq, such derivatives can lower the s-power by at most [(d + n)/2] units. 
We get non-vanishing contributions when n max = d + 2{k — r) > 0. They are proportional to the 
unperturbed coefficients 



Ak-r-n+[(d+n)/2];a(d+n)i ^-k-r-n~l+[(d+n)/2];a(d+n)+2 "' ^k~r-n;d+n ; (4-11) 



where n = 0, 1, • ■ ■ , ra max . In (|4.1ip A p - q is meant to vanish whenever p < 0. 

The classification applies to any unperturbed two-derivative operator Hq, in particular the 
spin-2 operator defined by gravity expanded around an arbitrary background. Finally, it can 
be easily generalized to operators Hq with a different maximal number of derivatives, to include 
fermions. 
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Non-minimal terms. Even if they are not multiplied by "small" parameters, non-minimal 
terms can be treated as perturbations, included in the n = term of Hi in (II. 2D . Indeed, each 
coefficient of the Schwinger-DeWitt expansion receives contributions from a finite number of non- 
minimal insertions and a finite number of commutators with Hq. Moreover, because non-minimal 
terms do not contain derivatives, their contributions are relatively easy to compute. For example 
to compute B 2 for 

H = □ + V, 

we can apply (|4.1ip with k = 2 and d = 0. We obtain non-vanishing contributions for n = 0, 1, 2, 
r = 1,2, namely a linear combination of V 2 Aq, VAq, VAi and VA^i. 



5 The case of gravity 

Expanding (II. ip around a background metric and choosing the harmonic gauge, the unperturbed 
spin-2 operator has the form 



H 0iJ a = a [ P 2^ P a ~ P 0vv ° ) + nonminimal terms, 

where 

p A' _l( S p' S <r' . §P ' S *' _ I a P'A p P'°' _ 1 Q PV 

are the projectors on the traceless and trace components, respectively. Define the bitensor 
H Q £ a (s;x,x') as the solution of 

i§;H /"\s;x,x>) + H 0p fH J"'(s-,x,x>) = 0, (5.1) 
with the boundary condition 

H o /°'(0;x,x') = ^f + ff s^Hx-x'). (5.2) 

Write the Schwinger-DeWitt expansion of H £ a (s;x,x') as 

H J°'(s;x,x') = "(^exp (^g^) Y.^) n A n ^\x, x>). 

The most general higher-derivative perturbation can be written as 

oo 

K^' = E V " V " tt ~"»23 M1 • ■ • (5.3) 

n=0 

where V^u pa ' W M ™ are completely symmetric tensors in the indices /ii • • • fi n , while the other 
indices satisfy obvious symmetry properties. The CBH approach described in this paper can be 
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applied with virtually no change. For example, in the case of a single insertion formula ([37 
generalizes to 

T(1) = \f ds J a i xJ^)[H 1 /jH /»\s- X , X ') 



and (|3,10p becomes 

oo [n/2] , , 



div 647r 2 e / y 4fcyM( n _ ^P^i ^n-afc k 



n=0 fc=0 
pa\fii---p n 



For renormalization purposes V^, can be taken to be traceless in [i\ ■ ■ • which amounts 

to exclude terms proportional to the field equations in (j5.3|) . Then formula (13. lip becomes 



1 r °° 

r (1) = — — / d 4 ^Vi„.; v ' v , r h ""'" ! 

n=0 



6 Conclusions 

In this paper we have studied improved Schwinger-DeWitt techniques for higher-derivative per- 
turbations of operator determinants and Green functions, to calculate counterterms and short- 
distance expansions of Feynman diagrams. In the perturbative regime the differential approach 
presents some difficulties, but it can be efficiently superseded by a systematic use of the CBH 
formula. We have classified the contributions that arise in this framework and outlined a number 
of simplification techniques. In some cases the calculational effort reduces considerably, in par- 
ticular when total derivatives can be neglected. The procedure is very general and applies also 
to quantum gravity treated with the background field method. 

Certain identities, such as (|3.9p . are new results, to our knowledge. They have been used 
to derive the closed formulas (|3.10p and (|3.1ip that relate the most general single-insertion per- 
turbed Schwinger-DeWitt coefficients to the unperturbed ones, up to total derivatives. Another 
new result is the three-derivative one- loop perturbed coefficient (|4.10p . When total derivatives 
are included and/or multiple insertions are considered, the list of contributing unperturbed co- 
efficients becomes considerably long. Nevertheless, we point out the simplicity of formulas (|4,7p 
and (|4.1U|) . compared with the involved intermediate expressions that lead to them. In particular, 
the inclusion of total derivatives in (|4.10p does not make the result much more complicated than 
(|3.8p . because several terms of (|3.8p are canceled by the total-derivative contributions. These 
facts suggest that there should exist more powerful and systematic simplification methods than 
the ones uncovered here. Hopefully the techniques of this paper can be extended and combined 
with the background field method to study two-loop and higher-order radiative corrections. 
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A Appendix: conventions 

The conventions used in this paper are those dubbed "SecondUp" (i.e. the default ones) in the 
package Ricci [T7] with metric signature (+, +, +, — ). Precisely, if is a vector, 

^p\vp Vp-.pv — RpovpV , Rpu — R iivpi R — Ftp,' 

We have performed our computations with two independent methods. The first method used a 
Mathematica package written by one of us (D.A.), the second method used the Ricci package. 
The first few Schwinger-DeWitt coefficients in the coincidence limit are 



^0 — 1, ^-0;p — 0, Aq-^v — -Rpv, A\ — — R, Ai-n — — — — R^. 
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